Anharmonic effects in MgB 2 ? 
A comparative inelastic X-ray scattering and Raman study 
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We study anharmonic effects in MgB2 by comparing fnelastic X-ray and Raman scattering to- 
gether with ab-initio calculations. Using high statistics and high q resolution measurements we 
show that the E2 g mode linewidth is independent of temperature along YA. We show, contrary to 
previous claims, that the Raman-peak energy decreases as a function of increasing temperature, a 
behaviour inconsistent with all the anharmonic ab-initio calculations of the E2 9 mode at Y available 
in literature. These findings and the excellent agreement between the X-ray measured and ab initio 
calculated phonon spectra suggest that anharmonicity is not the main mechanism determining the 
temperature behaviour of the Raman-peak energy. The Raman E2 9 peak position and linewidth can 
be explained by large dynamical effects in the phonon self-energy. In light of the present findings, 
the commonly accepted explanation of the reduced isotope effect in terms of anharmonic effects 
needs to be reconsidered. 

PACS numbers: 74.70.Ad, 74.25.Kc, 63.20.Dj, 63.20. Kr, 78.70.Ck, 71.15.Mb 
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I. INTRODUCTION 

The discovery of MgB 2 4 demonstrated that electron- 
phonon mediated superconductivity can generate 
critical-temperatures (T c = 39K) similar to those of 
perovskite oxides superconducting systems such as 
cupratesS. and doped BaBiOg^. Almost as surprising 
as the 39 K critical temperature is the reduced isotope 
effect^. The partial isotope coefficients a(X), defined 
as a(X) = — -rr^rf where Mx is the atomic mass of 
ion X, are a(B) =0.30(1) and a(Mg) = 0.02(1). Thus, 
the total isotope coefficient of a w 0.32 is substantially 
reduced from the canonical BCS value of cxbcs = 0.5. 

Deviations from olbcs value are usually attributed to 
strong electron-phonon coupling and to Coulomb repul- 
sion. However, theoretical investigations of the isotope 
coefficient using single- or double-gap Migdal-Eliashberg 
equations^ lead to unsatisfactory results since the iso- 
tope effect reduction is substantially underestimated. For 
this reason, it was proposed&i 8 -! 9 - that the T c in MgB 2 
is determined by the interplay of electron-phonon cou- 
pling and anharmonic effects. The latter can reduce T c 
through a hardening of the phonon frequencies. Refsi&£ 
showed than a hardening of « 25% of the E 2s phonon- 
frequency, along the high symmetry direction TA, would 
be necessary to explain the reduced isotope-coefficient. 

Anharmonic effects contribute to the broadening of a 
phonon and determine energy shift from its harmonic 
value. The anharmonic shift can be determined from the 
actual experimental phonon energy and from the har- 
monic energy, obtained from e.g. calculations. The in- 
trinsic linewidth of a phonon is, in general, given by two 
contributions due to anharmonicity and electron-phonon 



coupling. The anharmonic linewidth can, thus, be de- 
termined from the actual experimental linewidth and 
from the electron-phonon linewidth obtained e.g. by cal- 
culations. Finally, the temperature dependence of the 
phonon energy and linewidth is usually almost entirely 
due to the anharmonicity. The experimental determina- 
tion of the temperature dependence of the phonon energy 
and linewidth thus provide a further source of informa- 
tion. 

Raman spectra of MgB^Li^^J^l^il show a sin- 
gle very-broad peak corresponding to an excitation with 
E 2s symmetry at 75 meV, (i.e. 15% higher than the the- 
oretical harmonic frequency of 65 meV). Furthermore, 
the linewidth of the peak has a very strong temperature 
dependence. These two facts have been interpreted by 
several authors&^ i 12 ' 18 ! 19 as a signature of strong an- 
harmonicity. The presence of strong anharmonic-effects 
was further supported by ab-initio calculations based on 
the frozen-phonon model^&^ i 12 ' 15 : 18 : 19 . However, the re- 
sults obtained with this model are questionable, since the 
frozen-phonon approximation is an oversimplified model 
which does not include all the terms of the anharmonic 
perturbation theor y 20 : 21 : 22 : 23 : 24 . 

Ab-initio calculations of anharmonic effects includ- 
ing all the leading terms in anharmonic perturbation 
theory ( including scattering between different wave- 
vectors and different phonon modes ) 23 : 25 show that the 
phonon frequency shift of the E 2ff mode is only 5% of 
its harmonic phonon frequency at T. This is the re- 
sult of a partial cancellation between a large positive 
four-phonon scattering term (« +10%) and a negative 
three-phonon scattering term (« —5%). This is differ- 
ent from what is found by the more-approximate frozen- 
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phonon calculations, where the four-phonon scattering 
term is dominant^ 8 -! 9 - ' 12 i 15 ' 18 i 19 , This difference should be 
reflected in the temperature dependence of the Raman- 
peak energy, since the four-phonon scattering term gen- 
erates an increase in energy with temperature, while the 
sum of three- and four-phonons term should be slightly 
negativ e 23 ' 24 . 

Inelastic X-ray scattering (IXS) measurements do 
not support the hypothesis of strong anharmonic ef- 
fects. In fact, the phonon dispersion measured by 
is in good agreement with harmonic-frequency 
calculation !] 25 ' 27 and is consistent with a small anhar- 
monic shift along the TA direction. In particular, the 
IXS E2 9 phonon frequency near T is 65 meV, i.e. ~10 
meV smaller than the Raman peak. Ab-initio calcula- 
tions of the phonon linewidth are in very good agree- 
ment with IXS measurements^ 5 - and give an anharmonic 
contribution to the phonon linewidth which is negligi- 
ble with respect to the larger electron-phonon coupling 
contribution. Moreover, calculations^ 5 - reveal a very 
weak temperature-dependence of the phonon-linewidth 
between 50 K and 300 K. Finally, the authors of Ref, 2e 
claim that the phonon-linewidth of the F,2 g mode is in- 
dependent from temperature between 40 K and 300 K 
along the TA. If confirmed, this last finding would es- 
tablish that anharmonic effects are weak. However, the 
moderate resolution and the fact that in Refi^ the F,2 g 
mode is seen only as a shoulder of the nearby Ei u mode 
do not allow for a definitive conclusion. 

In this work we present a high-resolution inelastic X- 
ray scattering and Raman study, in order to settle the 
debate on the presence of important anharmonic effects 
in MgB2. The two techniques are complementary since 
inelastic X-ray scattering probes the phonon dispersion 
throughout the entire Brillouin zone. However, at the 
r point the elastic (Bragg) scattering usually dominates, 
therefore masking the phonon contributions, as it is the 
case in MgB2- On the contrary, Raman experiments 
probe excitations at very small exchanged momentum, 
q ~ 0. We show with high statistics and high q-rcsolution 
IXS measurements that the E 2g mode linewidth is inde- 
pendent of temperature along T-A, and furthermore de- 
termine the dispersion along T-M with the best resolution 
to date. We investigate the behaviour of the Raman-peak 
energy and linewidth as a function of temperature. All 
our results are consistent with the presence of small an- 
harmonic effects in MgB2. We, finally, show that, on 
the basis of recent theoretical result s 24 ' 28 ' 29 , there is no 
contradiction between the presence of a Raman-peak at 
75 meV and the X-ray E 2g phonon frequency at 65 meV 
near T. 

The paper is structured as follows. In sec. [Til we de- 
scribe the methodological details involved in inelastic X- 
ray scattering, Raman spectroscopy and ab-initio simu- 
lations. In sec. Mil we present our main results and in 
sec. IIVI we comment on the results and discuss the main 
implications for the interpretation of Raman data. 



II. METHODS 
A. Experiments 

1. Sample 

The crystal used in our experiment was grown at a 
pressure of 30-35 kbar. A mixture of Mg and B was put 
into a BN container in a cubic anvil device. The temper- 
ature was increased during one hour up to 1700— 1800°C, 
kept stable for 1-3 hours, and finally decreased during 1-2 
hours. As a result plate-like MgB2 crystals were formed 
of which we used a sample of about 1.00 x 0.45 x 0.08 mm 3 
with a measured in-plane mosaic of 0.011°. 



2. Inelastic X-ray Scattering 

The experiment was carried out at the Inelastic X- 
ray Scattering Beamline II (ID28) at the European 
Synchrotron Radiation Facility in Grenoble (France). 
The x-ray beam from the undulators was monochro- 
matized by a cryogenically cooled silicon (111) double 
crystal monochromator, followed by a back-scattering 
monochromator, operating at a Bragg angle of 89.98 de- 
grees, and utilizing the Si (9 9 9) reflection order—. The 
back-scattered beam was then focused by a platinum- 
coated toroidal mirror, which provided a focal spot at 
the sample position of 0.270 (horizontal) and 0.080 mm 2 
FWHM at the sample position. The scattered pho- 
tons were energy-analyzed by a 6.5m Rowland circle 
five-crystal-spectromete r 31 ' 32 . The analysers operate at 
the same reflection order as the monochromator, corre- 
sponding to a wavelength A = 0.696782 A -1 (photon en- 
ergy: 17794 eV). The overall instrumental function has 
a pseudo-Voigt line-shape with a width AE of 3 meV 
(FWHM) as determined from the elastic scattering of a 
Plexiglas sample (at a Q-transfer of 10 nm^ 1 , and T — 
10 K). The energy-analyzed photons are detected by a 
Peltier-cooled silicon diode detector which has an intrin- 
sic energy resolution of 400 eV. The dark counts due to 
electronic and environmental noise amounts to about 0.2 
counts/minute. Further components of the spectrometer 
are an entrance pinhole, slits in front of the analysers in 
order to define the momentum transfer resolution and a 
detector pinhole for aberrant ray suppression. 

The analyser crystals are temperature stabilized at 
295.65 K with a typical stability of 1 mK/24h. The 
momentum transfer Q — 2 ki sin(|), where ki is the 
incident photon wave vector and 8, the scattering angle, 
is selected by rotating the spectrometer around a vertical 
axis passing through the scattering sample in the hori- 
zontal plane. Since there are five independent analyser 
systems, spectra at five different momentum transfers can 
be recorded simultaneously. Their separation is energy 
dependent and for the Si(9 9 9) reflection it amounts to 
2.43 nm -1 . The energy scans are performed by varying 
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the lattice spacing of the monochromator temperature 
while the analyser temperature is kept fixed. Conversion 
from the temperature scale to the energy scale is accom- 
plished by using the known thermal expansion coefficient 
of silicon: a(T) = a + (T-T ) with a = 2.581 x 1CT 6 
1/K and 13 = 0.016 x 10~ 6 1/K 2 , T = 295.65 YM. 

The validity of this conversion has been checked by 
comparing the measured diamond dispersion curve for 
longitudinal acoustic and optic phonons with well estab- 
lished inelastic neutron scattering results. 

The energy scans were fitted using a sum of pseudo- 
Voigt: 



r/2 



(e-e ) 2 +r2/4 
and Lorcntzian functions: 



ryexp 



(ln2)(e-e ) 1 

r 2 /4 



r/2 



(e-e ) 2 +r 2 /4 



(1) 



(2) 



functions, where e = huj = hv is the energy. Pseudo- 
Voigt functions were used to fit elastic and resolution- 
limited inelastic contributions, with T and r\ parame- 
ters adapted to match the instrumental function, while a 
Lorentzian line-shape was used to fit the E 2g mode, with 
a width r fitted to the data. Alternatively, a convolution 
with the experimental resolution function has been used. 
In this case, the response function was composed of a sum 
of delta functions (for elastic and resolution-limited in- 
elastic contributions) plus a Lorentzian for the E2 9 mode. 
The two data analysis procedures gave consistent results. 

The fitting algorithm used a ^-minimisation 
routine 34 , with the condition that the detailed balance 
between Stokes and anti-Stokes excitations is fulfilled. A 
constant background, coming essentially from electronic 
noise, was added. 

The MgB2 crystal was glued on a copper sample 
holder, which was mounted inside a vacuum chamber 
on the cold finger of a closed-loop helium cryostat (Cry- 
omech ST15). The temperature was measured by a sil- 
icon diode attached to the sample holder. An indepen- 
dent check of the sample temperature was provided by 
the Stokes- Anti-stokes intensity ratio. 



3. Raman measurements 

Raman spectra were excited with the 488 nm line of 
an Ar laser focused with an optical microscope (long dis- 
tance lOOx objective). The scattered light was analyzed 
by a single-grating Labram spectrometer equipped with 
a CCD. Temperature dependent measurements were per- 
formed in an Oxford flow cryostat with a sapphire win- 
dow to preserve the polarization of the light. We used 
a 600 1 /nm grating to record a wide spectral window 



(5000 cm -1 or 620 meV). The loss in resolution is irrel- 
evant for MgB 2 , because its Raman features are at least 
several hundreds of wave numbers in width. Some ex- 
periments were also performed with 514 nm and a triple- 
grating spectrometer to compare with published studies. 
We report the as-measured data normalized to the exci- 
tation energy and accumulation time. Unless noted oth- 
erwise, the spectra were not scaled or shifted along the y 
axis. 



To analyze the symmetry of the scattered light we 
recorded Raman spectra in various backscattering ge- 
ometries kiie-i, e s )k s , where fc, and k s indicate the di- 
rection of the incoming (i) and scattered (s) light and ej 
and e s the polarization (Porto's notation) 35 . For fej||fc s ||c 
we used backscattering normal to the surface; in the 
ki\\k s \\a,b configurations the laser was focused onto the 
side of a MgB2 crystal. The polarization of the incoming 
light was rotated by a Fresnel rhomb. The polarization of 
the scattered light was chosen by a A/2 wave plate com- 
bined with an analyzer. We verified the polarization- 
dependent measurements by recording spectra on the 
(111) and (001) surface of silicon and on single crystals of 
graphite, which has the same point group as MgB2. We 
also ensured the internal consistency of the Raman inten- 
sities; e.g., the spectra recorded in c(a,a)c and b(a,a)b 
were identical in shape and intensity as required by the 
Raman selection rules. 



The frequency of the ~ 600 cm -1 (74 meV) peak of 
MgB2 reported in the literature varies between 74 and 
77 meV at room temperature! 10 ' 13 i 14 ' 16 i 36 Part of the dis- 
crepancy might be due to a varying sample quality. How- 
ever, a major contribution comes from the uncertainty re- 
sulting from the functional form used for the background 
and the peak. We tested several fitting routines and 
found that the uncertainty for the position and width 
of a Lorentzian is around 1.5% or 1.2 meV. On the other 
hand, when using a Fano line-shape the bare phonon fre- 
quency can be as low as 71 meV or as high as 79 meV 
(300 K) depending on the sign of the asymmetry param- 
eter. The peak positions and line widths reported in this 
paper were obtained by assuming the 600 cm -1 peak to 
be Lorentzian and the continuous background in the Ei g 
spectra to follow the same frequency dependence as the 
ctxx(Aig) background signal (see Sect. IIIIBTl for details). 
Other fit functions shift the position and width. We also 
observed a change in the phonon frequency when reduc- 
ing the energy range in a fit. For our analysis we used 
the measured spectra between 100 and 2000 cm -1 (~ 10 
and 250 meV). Smaller windows yield a strong scattering 
in the peak positions. We stress, however, that varying 
fit routines resulted in a constant offset of the tempera- 
ture dependence or stronger scattering in the data. The 
functional dependence of the Raman peak position and 
width on temperature, however, remained unaffected. 
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B. Calculations 



Electronic structure calculations were performed using 
density functional theory (DFT) and phonon frequen- 
cies were calculated in the framework of linear response 
theory^ using the PWSCF/espresso package^. We used 
the generalized gradient approximation^? and norm con- 
serving pseudo-potentials^. For Mg we used non-linear 
core corrections^ and we treated the 2s and 2p levels as 
core states. The wave-functions were expanded in plane 
waves using a 35 Ry cutoff. The calculations were per- 
formed with the experimental crystal structure, namely 
a = 3.083 A and c/a = 1.142. For the phonon frequen- 
cies calculation, we used a 16 x 16 x 16 Monkhorst-Pack 
grid for the electronic Brillouin Zone integration and first 
order Hermite-Gaussian smearing^ of 0.025 Ry. 

The electron-phonon coupling contribution to the 
linewidth (full width half maximum ) at momentum q 
for the v phonon mode, with frequency cj q „ can be writ- 
ten as^: 



7q^ 



£ 



l3kn,k+gm| 2 <K £ kn)<K£k+gm) 



(3) 



where the k-sum is extended over the Brillouin Zone , 
Nk is the number of fc-points in the sum, and £kn are 
the electronic energies, measured with respect to the 
Fermi level, of the n band at point k. The matrix 
element is 5k„,k+qm = (kra|<JV/<Je qi/ |k + qm)/ v /2w q ^, 
where SV/Seq^ is the derivative of the Kohn-Sham po- 
tential V with respect to the phonon coordinate e q „ nor- 
malized over the unit cell. 

The electron energies £kn and the g£ n k+qm matrix ele- 
ment in eq. [3]are obtained fully ab initio. The actual sum 
over the Brillouin zone is calculated on a Nk = 27 x 27 x 27 
uniform k-point mesh shifted by a random vector with 
respect to the origin. 

The x-ray structure factor for one-phonon scattering 
at temperature T, frequency ui and reciprocal-space vec- 
tor q + G = Q (with q belonging to the first Brillouin 
Zone and G a reciprocal-lattice vector, and Q the total 
momentum transfer) is obtained from 



S(Q, w) = £ GM, G)F V (w, T, G) 



(4) 



with 



G,(q,G) 



F„(w,T,G) 



£/a(Q) — =[(Q). e «Je 



jQRo 



5{u + 7?W qi /) (5) 



where a labels the atoms in the unit cell, v the mode 
and / a (Q) is the atomic form factor calculated as in—. 
Moreover, M a , R Q and W a are the atomic mass, posi- 
tion in the unit cell and the Debye- Waller factor, respec- 
tively. Finally, (n ql ,) is the phonon occupation number 
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FIG. 1: (Color online) Positions of the analyzers 
squares) for the measurements along T-A (see Fig. [6]), 
near F (see Fig. 0. r=G(HKL)=(2,l,0). 



(red 
and 



and e^/x/Mf, is the atomic displacement of the atom a 
corresponding to the phonon eigenvector e q „. The value 
i] = — 1 (rj = 1) identifies Stokes (anti-Stokes) processes. 

In what follows we neglect the Debye- Waller factor and 
we smear the 6— functions with the convolution between 
the phonon intrinsic lineshape calculated ab initio and 
the instrumental function (3 meV FWHM). For the cal- 
culation of the x-ray structure factor we used the atomic 
factors for neutral Mg and B fromii and the phonon pat- 
terns calculated ab initio. 



III. 



RESULTS 



A. Inelastic X-ray Scattering results 

The reciprocal points in the following section are ex- 
pressed in terms of hexagonal parameters a = b = 3.0828, 
A c = 3.5186 A, a = (3 = 90° and 7 = 120°. In these 
coordinates, all the measurements were recorded in the 
Brillouin Zone centered around (2,1,0). 

In Fig. [1] and [2] we show the relative position of the 
analyzers for the measurements reported below. 



1. Dispersion along Y-M 

The experimental results shown in figure [3] were ob- 
tained exploiting the signal from all five analyzers. The 



5 




FIG. 2: (Color online) Positions of the analyzers (red 
squares) for measurements along T-M (see Fig. [IJ. 
r=G(HKL)=((2,l,0). 
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FIG. 3: (Color online) Left: IXS data (circle) and fit 
(lines) along F-M for T=46 K. Right: IXS data (circle) and 
calculations (lines) along F-M for T=46 K. The legend box 
indicate the f components of the q reduced vector (see text). 



relative Q displacement between the analyzers is fixed by 
the spectrometer geometry, and for longitudinal modes 
such a configuration can be close to a high symmetry 
direction. This opportunity has been as well exploited 
in previous works on MgB ^ 25 ' 26 . Here the difference is 
that, using a different incoming wavelength of 0.696782 
A -1 , we have higher energy resolution and larger steps 
in reciprocal space, corresponding to the different fixed 
angular displacement of the analyzers. For a position 
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FIG. 4: (Color online) Bottom panel: MgB2 phonon disper- 
sion. IXS data (circles) frorn^ (blue open) and the present 
work (black full), compared to calculations (red lines). Top 
panel: E2 9 mode width from fit of the IXS data of the present 
work (black diamonds) and Ref: 25 (blue open with errorbars) 
and calculations (red points with continuous line connected). 
Stars corresponds to the high energy and open circles to the 
low energy E2 9 modes when not degenerate. 



close to the T-M line, this implies Q points G + q = 
(2, 1, 0) + (£, v, C) where v and £ are small compared to 
£. The configuration chosen for this experiment was as 
follows: 



analyzer 5: 
0.104 r-M; 

analyzer 4: 
0.106 T-M ; 

analyzer 3: 
0.307 T-M; 

analyzer 2: 
.509 T-M; 

analyzer 1: 
0.700 T-M; 



: (2.052 1.000 0.0), |f | « 0.052 or 
: (1.948 1.002 0.0), |f | w 0.053 or 
(1.847 1.001 0.0), |f | w 0.1535 or 
(1.744 0.997 0.0), |f| w 0.2545 or 
: (1.645 0.990 0.0), |f | « 0.350 or 



The positions in the hexagonal plane are shown in Fig. 
[5] Note that analyzers 4 and 5 are placed symmetri- 
cally around the chosen zone center. These positions are 
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equivalent for phonon frequencies but not for dynamical 
structure factors. Contribution along < 0, v,0 > are less 
than 0.005 r.l.u. for analyzer 5 to 2 and 0.01 r.l.u. for 
analyzer 1. 

We measured the IXS spectra using a slit opening in 
front of the analyzers of 60 mm in both horizontal and 
vertical direction, amounting to a Q resolution of ±0.416 
nm -1 for analyzers 1, 2 and 3, while analyzers 4 and 5 
had 60 x 20 mm 2 vertical x horizontal opening. This 
corresponds to a main contribution of ±0.0022 r.l.u. in 
the ( , v , ) direction and to ±0.008 r.l.u. in the 
perpendicular ( , , C ) ( or T-A) direction. 

In Fig. [31 left panel, we show the IXS spectra obtained 
at ss 46 K for this configuration, together with the fit used 
in order to extract mode frequencies as well as linewidth 
values for the E 2ff mode. A recent paper— suggests that 
the MgB2 IXS spectra show additional structures due 
to two-phonon contributions. Although we do not have 
analyzed our data in this way, such contributions may 
actually explain small deviations of our fit results. In the 
right panel we show the same data compared with ab- 
initio calculated spectra, which correctly reproduce the 
measurements. 

The same analysis is done for points along T-A (see Sec. 
IIII A3j) . These are grouped together with previous re- 
sults at lower resolution 2 ^ in order to experimentally de- 
termine the phonon dispersion and the lincwidths along 
all the high symmetry directions. We can thus compare 
the experimental phonon dispersion with the theoretical 
calculation along the high symmetry lines, as shown in 
Fig. 0] bottom. In the top panel, we show the experi- 
mental line-width T, for the E2 S mode, de-convoluted by 
the instrumental broadening, in comparison with the cal- 
culated linewidth. The agreement of the calculated fre- 
quencies with experiment is remarkable. The same is ob- 
served for the agreement between the experimental line- 
width, de-convoluted by the instrumental function, and 
the calculated electron-phonon-coupling contribution to 
the intrinsic E2 9 linewidth. Note that data from Ref. 25 
are obtained at room temperature and with 6.1 meV in- 
strumental function FWHM. 



2. High resolution line-width Ei g measurements close to F 

Fig. [5] shows the room temperature MgB2 IXS spec- 
tra for a Q-value corresponding to 0.08 xT-A, together 
with the data, we show calculation convoluted with the 
instrumental resolution for all modes, including the E 2s 
one. The calculation was scaled as in Fig. [51 

We measured the IXS spectra using a slit opening 
in front of the analyzers of 20 mm in both horizontal 
and vertical direction, equivalent to a Q resolution of 
±0.139 nm . This corresponds to a main contribution 
of ±0.0022 reduced lattice units (r.l.u.) in the ( , v , 
) direction and to ±0.0026 r.l.u. in the perpendicular ( 
, , C ) ( or T-M) direction. 




60 80 
Energy (meV) 
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FIG. 5: (Color online) Room temperature IXS data at 
0.08 xF- A. The line is an ab-initio simulation convoluted with 
the energy resolution, without electron-phonon coupling 
broadening effect on the E29 mode. 



Temperature- dependence of the IXS spectra 



In figure [S] we compare IXS E2 S spectra measurements 
and calculations performed at two different sample tem- 
peratures (~ 50 K and 300 K) for Q=(HKL) correspond- 
ing to a phonon propagation vector of (0, 0, 0.3) (top 
panel) (i.e. at 60% of the T-A line) and the zone bound- 
ary A point (0,0, 0.5) (bottom panel). 

To induce a maximum effect from phonon-phonon scat- 
tering (anharmonicity) 23 ' 25 , we choose the largest pos- 
sible temperature difference in the normal-state region, 
i.e. from room temperature to T above ~ 40 K. This 
was done to avoid effects due to the superconducting gap. 
More precisely, the measured temperature by the probe 
on the sample holder was stable at 46 K and 298 K re- 
spectively. 

We measured the IXS spectra using a slit opening 
in front of the analyzers of 60 mm in both horizontal 
and vertical direction, equivalent to a Q resolution of 
±0.416 nm -1 . This corresponds to a main contribution 
of ±0.0065 r.l.u. along the ( , v , ) direction and to 
±0.008 r.l.u. along the perpendicular ( , , C ) direction 
(i.e. along T-A). The theoretical spectra are arbitrarily 
normalized to reproduce the height of the highest vibra- 
tional peak (Ei u , not shown in Fig. [5]) 
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FIG. 6: (Color online) Left: IXS data at T=300 K (red full 
points, with errorbars) and T=46 K (blue empty points). 
Right: IXS data at T=46 K (blue empty points) and ab-initio 
simulation (lines) Measured spectra and calculation spectra 
at 0.6 xF- A (top) and at the A point (bottom). 



FABLE I: Raman selection rules for the D^h point group of 
MgB2 . a and b are vectors within the hexagonal plane forming 
an angle of 90° . Since the hexagonal lattice is isotropic in the 
plane a and b can be chosen arbitrarily (with respect to the 
lattice vectors), c is normal to the plane. Fhe an denote the 
element ii of the corresponding Raman tensor. 
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Table U gives the selection rules for the Raman active 
representations of MgB2 {Dqu point group). Note that 
the intensity of E^ g is independent of polarization as long 
as the light incidents and scatters normal to the ab plane. 
In Fig. [7] we present the four independent components of 
the Raman scattered light from the MgB 2 crystal. The 
~ 75 meV peak is indeed only found in the Ei g scattering 



1500 
1000 
500 


J2 1500 

'E 

_q 1 000 

500 

'in 

c 

CD 

£ 

1500 
1000 
500 



1 1 1 l: 
- \a> 

j 1X5 j 

mm: 


V 1 1 1 | 1 1 1 1 | 1 1 1 1 

: ioo so 

120^^-——,-^ 60 

■ \ iso l — ; m r^T^^V: — * — Jo - 




i Raman £, g [c{a,b)c] '- 




t i i i | i i i i | i i i i 

E uj [b(c,a)b] ; 


E [b(a,c)b] ; 

{•'• \ > 'i > 'i 


" (c) 

'-_ (C5T 


^-^A uj [b(c,c)b] 


u \ 


A^[c{a,a)c-c{a,b)ch 

■.. \ .... \ .... \ ... r 



50 100 150 200 

Raman shift (meV) 

FIG. 7: (Color online) Raman selection rules and symmetry 
of the scattered light in MgB2 (room temperature), (a) Ei g 
component with the peak at 75 meV (608 cm -1 ). The back- 
ground measured under crossed polarization is also of Ei g 
symmetry in contrast to parallel polarization, where it has 
E2 g and Ai g [see (c)] contribution. The inset shows the con- 
stant integrated intensity of the 75 meV peak when rotating 
the in-plane polarization, (b) Ei g symmetry obtained in two 
scattering configurations, (c) a zz and a mx ~ ctyy A\ g com- 
ponent of the scattering intensity. The latter was obtained 
by subtracting the c(a, a)c and c(a, b)c Raman spectra. The 
two dashed lines mark the phonon frequency of the E 2g peak 
measured with Raman and inelastic X-ray scattering. 



configuration , 10 ' 13 Fig. Ufa), and its intensity is indepen- 
dent of the in-plane orientation of the crystal (see inset). 
The weak feature at 75 meV in the E\ g b(a, c)b configu- 
ration is within experimental error. The two dotted lines 
mark the En g peak frequency as obtained by Raman (75 
meV) and inelastic X-ray scattering (w 69 meV) at 0.1 T- 
M and ab-initio calculations at T (70 meV)^ 3 -. As noted 
earlier, the two measurements give phonon frequencies 
differing by at least ~ 6 meV. 

Our Raman spectra, consistently with the 
literatureiS'ii^ii 3 -^^^'^, are smooth. They do 
not show any structure that can be interpreted as 
a multiple-phonon peak, as proposed in Ref. [45l In 
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FIG. 8: (Color online) E2 9 Raman spectra (raw data) of MgEb 
between 5 K and room temperature in c(a, h)c scattering con- 
figuration. An offset has been added for clarity (see the zero 
line at the cut-off of the notch filter). 



general, two-phonon Raman spectra have more than 
one symmetry component^. However the A\ g and E\ g 
spectral contributions in Fig. do not show any sharp 
feature. 

The polarization dependence of the continuous back- 
ground for a xz (E lg ), Fig.Efb), and a zz (A lg ), Fig.Efc), 
agree with Quilty et aZ.— . The scattering intensity in- 
creases linearly with the frequency of the excitation. We 
extracted the a xx (Ai g ) background component from the 
measured spectra as shown in Fig. 0(c). In contrast to 
the frequency dependence of the zz and xz component of 
the Raman tensor, a xx has a maximum at ~ 60 meV and 
decreases in intensity towards smaller frequencies. The 
a xy (E 2g ) background component appears to be similar 
to a xx (Ai g ); the overall intensity is almost identical and 
the decrease in background intensity towards small fre- 
quencies is present as well as can be seen in Fig. 0(a). 
This points towards a common origin of the electronic 
Ei g and A\ g scattering, e.g., such a behavior might arise 
from a coupling of electrons that are correlated with the 
E\ g representation of the MgB 2 point group. 

The precise origin of the electronic background in the 
MgB-2 Raman spectra is not known. It is remarkable that 
the signal is strong and approximately constant over a 
wide energy range (at least up to 600 meV the highest 
energy in our measurements). 



ence of the Raman spectra 



2. Temperature- 



Figure [8] shows Raman spectra measured between 5 K 
(bottom) and room temperature (top). The 75 meV 
(600 cm -1 ) peak decreases in frequency and broadens at 
higher temperatures. The frequency shift is more pro- 
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FIG. 9: (Color online) Position (a) and full width at half 
maximum (b) for the 75 meV (600 cm -1 ) Raman peak. The 
spectra were fitted by a Lorentzian plus a background. The 
frequency dependence of the background is assumed to be the 
same as the a xx (Ai g ) component, see Fig. [re- 



nounced than reported by Quilty et al^ and disagrees 
with the finding by Martinho et a/.—, who reported a 
constant peak frequency between 10 and 300 K. The shift 
in the peak maximum agrees with Rafailov et al*£, al- 
though this is somewhat difficult to judge given the nar- 
row Raman frequency window in the reference. Surpris- 
ingly, we find the electronic background to be indepen- 
dent of temperature. The intensity is constant (the spec- 
tra in Fig. [8] are raw data) and there is no increase in the 
scattering intensity at low frequency for low r 10 i 36 . 

The frequency of the ~ 75meV (600 cm -1 ) Raman 
peak decreases from 81 meV to 75 meV between 5K and 
room temperature as we show in Fig. [Ha). This is ac- 
companied by an increase in the line width by 45% or 
from 37 to 54 meV (Fig.^b)). 



IV. 



DISCUSSION 



The difference between the 75 meV Raman peak and 
the calculated harmonic F,2 g phonon frequency (65 meV) 
has been attributed to the presence of anharmonic effects 
which are supposed to be large& along the entire TA direc- 
tion. Here, we have determined using high resolution IXS 
the phonon dispersion and the phonon linewidth along 
three high symmetry directions in the Brillouin zone to 
verify whether anharmonic effects are relevant or not. 
As already shown in Refsj 25 ' 26 , the IXS phonon disper- 
sion is in good agreement with the calculated harmonic 
phonon dispersion. To further determine the magnitude 
of anharmonic effects, in the present work, we have also 
measured the phonon frequency shift and linewidth as a 
function of temperature using Raman spectroscopy and 
inelastic X-ray scattering. 

Our Raman spectra show a well defined peak at « 
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£ = 0.0 




= 0.3 


£ = 0.5 




Raman 


Theory 
(Anharm.) 


IXS 


Theory 
(EPC) 


IXS 


Theory 
(EPC) 


Theory 
(Anharm.) 


T=300K 


53 ± 3 


1.21 


23 ± 2 


23.0 


26 ± 2 


20.35 


2.13 


T=50K 


39 ± 2 


0.16 


20 ± 2 


23.0 


24 ± 2 


20.35 


0.10 



TABLE II: Measured (Raman and IXS) and calculated phonon-linewidth (FWHM) of the E2 S phonon mode at the (0, 0, £) 
wave-vector in meV. Both the electron-phonon coupling (EPC) and the anharmonic contributions (Anharm.) to the phonon 
linewidth are reported. Computational details are explained in Refs . 23 ' 25 . Note that the anharmonic contribution has been 
calculated only at selected high symmetry point: F (£ = 0.0) and A (£ = 0.5). There is no EPC contribution at T^-. 



75 meV having E 2g symmetry. The energy of the 75 
meV peak shifts downward and broadens by increas- 
ing the temperature. The frequency shift is more pro- 
nounced than reported by Quilty et alJ^ and disagrees 
with the finding by Martinho et al^, who reported 
a constant peak frequency between 10 and 300 K. The 
decrease in energy of the Raman-peak on temperature 
increase is not consistent with the claim that anhar- 
monic four-phonon scattering is the dominant anhar- 
monic contribution 6 -* 8 -! 9 - ! 12 ! 18 ' 19 . Indeed, since the four- 
phonon scattering term is positive at T— & & 12 ' 18 i 19 ' 23 , 
the corresponding anharmonic shift should increase with 
temperature (see e.g Eq. 4 in Refi^ 3 - ). However, even 
including a negative third order term, state-of-the-art 
calculations 2 - report a substantially smaller decrease of 
the F,2 g phonon frequency than what is presently mea- 
sured. This suggests that anharmonicity is not the main 
mechanism determining the temperature behavior of the 
Raman-peak energy and that the Raman-peak is not due 
to a bare phonon-vibration. We also remark that the vari- 
ation of the Raman shift with the temperature (Fig. [5]) 
cannot be explained by a two-phonon contribution as pro- 
posed by Refl45l (see e.g. Ref. [H). 

The large increase of the Raman linewidth between 
40 K and 300 K has no equivalent in the IXS spec- 
tra. Using high-resolution IXS we have shown that the 
phonon linewidth along the TA direction is essentially 
temperature-independent in the 46-300 K range. This 
was already claimed in a previous paper—. However 
in Reb2£ the E 2ff peak was only detected as a shoul- 
der of the nearby E ltl mode, while in the present work 
we have a very well resolved peak for the E2 9 mode in 
the region from 0.5rA to A. As shown in Table HT1 and 
Fig. [5] (top panel), ab-initio simulations are successful 
in predicting the phonon linewidth (including electron- 
phonon and anharmonic damping) of MgB2 along all the 
high-symmetry directions considered. This very good 
agreement confirms that, in the region from 0.5TA to A, 
phonon-damping is almost exclusively due to electron- 
phonon coupling and not to anharmonic effect s 23 ' 25 . 

We finally want to stress that there is no contradic- 
tion between the Raman-peak at 75 meV at T and the 
X-ray E2 ff phonon frequency at 65 meV near T and ab- 
initio harmonic calculations (giving 65 meV for the E2 S , 
r phonon). The key to understand these differences is 
the presence of relevant dynamical and electronic effects 



in the phonon self-energy, as shown in Ref i ' i ' . So 
far, ab-initio phonon calculations in MgB2 have been per- 
formed within the adiabatic approximation. The adia- 
batic approximation assumes that the phonon is a static 
perturbation, but, actually, a phonon is dynamic pertur- 
bation, oscillating at the frequency uo. Within the adi- 
abatic approximation, the phonon self-energy does not 
depend on u>, while, in general, the phonon self-energy 
depends on lu. The adiabatic approximation is valid if 
the interatomic force constants are istantaneous, that is 
the force on a given atom at a given time depends on the 
position of the other atoms at the same time^l. In real- 
ity, however, the force on a given atom at a give time t 
depends on the position of the other atoms at a previous 
time t— . This dynamical effect can be treated within 
time-dependent perturbation theory and can be relevant 
for zone-center optical phonons. 

As shown in Refs i 24 ' 28 , the F<2 g phonon frequency 
varies strongly for q near T due to dynamical effects in 
the real-part of the phonon self-energy. As an example, 
along the TM direction, the E2 S phonon frequency has a 
smooth behavior for q > 0.1TM and increases abruptly 
for q < 0.05TM (see fig. 3 in Refs^i or fig. 3(a) of 
Ref.—). The E2 9 X-ray dynamical structure factor is 
strong enough for accurate statistics only for q > 0.1TM, 
that is IXS are probing the smooth region. On the con- 
trary, Raman spectroscopy probes excitations at much 
smaller q (q < 0.05TM), i.e. in the region where the 
phonon frequency is higher. Concluding, the w 10 meV 
difference between the Raman and X-ray peaks is ex- 
pected from the theoretical findings of Rcfs. 24 ' 28 and 
is due to the pathological behavior of the phonon self- 
energy at r. 

The disagreement between Raman data and the har- 
monic ab-initio calculations, published so far, is due to 
the fact that these calculations have been performed 
within the adiabatic approximation. When dynamical 
effects are neglected, the mentioned pathological behav- 
ior is absen t 24 ' 46 and the phonon dispersion is smooth 
also in the region near T. Therefore, the calculated F<2g 
frequency at T is smaller than the actual frequency, which 
should be computed including dynamical effects. These 
dynamical effects have no influence outside a small region 
around T and, for this reason, the ab-initio harmonic cal- 
culations done outside this region are correct and repro- 
duce the phonon dispersion. 
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Finally, in Ref^ it was shown that the Allen formula^ 
gives a zero linewidth for the E 2s mode at T. As shown 
in fig [5] this is not the case. To explain the presence 
of a non-zero linewidth at T, it is necessary to include 
electron self-energy effects in the phonon self-energy as 
pointed out by Cappelluti in Ref. 29 . According to Rcf. 29 
these effects are responsible for the unexplained softening 
of the Raman-peak energy and increase of the Raman- 
peak linewidth by increasing the temperature. However, 
the model of Ref3?> is simplified, therefore giving only a 
qualitative trend. More theoretical work is required to 
achieve quantitative agreement with experiments. 

V. CONCLUSIONS 

We have presented a high-resolution inelastic X-ray 
scattering and Raman study, in order to settle the debate 
on the presence or not of important anharmonic effects in 
MgB2. First, we have shown with high statistics and high 
q resolution measurements that the E2 S mode linewidth 
is independent of temperature along TA, ruling out a ma- 
jor contribution of anharmonicity in the lattice dynamics 
of MgB2. Moreover we measured the dispersion along T- 
M with the best resolution to date. The E2 9 mode shows, 
along this line, large variation of both phonon frequencies 
and line-width, as previously observed along A-L 2 ^ and 
T-M 2 ^ directions. Second, we have investigated the be- 
havior of the Raman-peak energy and linewidth as a func- 
tion of temperature. The Raman-peak energy decreases 
as a function of increasing temperature. This behavior is 
not reproduced by anharmonic ab-initio calculations of 
the E 2ff T mode (calculations done both with the over- 



simplified frozen-phonon approach£&£ i 12 i 18 i 19 i 23 or with 
state-of-the-art perturbation theory 2 ^). This finding sug- 
gests that anharmonicity is not the main mechanism de- 
termining the temperature behavior of the Raman-peak 
energy and that the Raman-peak is not due to a bare 
phonon- vibration. Finally, we have shown that, on the 
basis of recent theoretical result o 24 ' 28 i 29 , there is no con- 
tradiction between the presence of a 75 meV Raman-peak 
at r and the X-ray E2 fl phonon frequency at 65 meV near 

r. 

In conclusion, the present results indicate that anhar- 
monicity plays a marginal role in MgB2. As a conse- 
quence, the explanation of the reduced isotope effect*^ 
- one of the most important unresolved issues in the 
physics of MgB2 - needs to be reconsidered. At present, 
the explanation of the reduced isotope effect is one of the 
most important unresolved issues in the physics of MgB2 . 

Acknowledgments 

We wish to acknowledge the ESRF for the support of 
experiment HS-2598. The Raman measurements were 
performed at the Technische Universitat Berlin. We 
thank C. Thomsen for kind hospitality in his lab, the use 
of the Raman equipment, and helpful discussions. We 
are grateful to S. Bahrs and D. Heinrich for help with 
the temperature-dependent Raman measurements. We 
thank A. Ferrari for discussion on the Raman measure- 
ments. Calculations were performed at the IDRIS super- 
computing center (project 071202). This work was sup- 
ported by the Swiss National Science Foundation through 
NCCR MaNEP. 



1 J. Nagamatsu, N. N. end T. Muranaka, Y. Zenitani, and 
J. Akimitsu, Nature 410, 63 (2001). 

2 J. G. Bednorz and K. A. Miiller, Z. Phys. B 64, 189 (1986). 

3 R. J. Cava, B. Batlogg, J. J. Krajewski, R. Farrow, L. W. 
Rupp, A. E. White, K. Short, and W. F. P. T. Kometani, 
Nature 332, 814 (1988). 

4 D. G. Hinks, H. Claus, and J. D. Jorgensen, Na- 
ture(London) 411, 457 (2001). 

5 S. L. Bud'ko, G. Lapertot, C. Petrovic, C. E. Cunningham, 
N. Anderson, and P. C. Canfield, Phys. Rev. Lett 86, 1877 
(2001). 

6 H. J. Choi, D. Roundy, H. Sun, M. C. ML, and S. Louie, 
Nature (London) 418, 758 (2002). 

7 H. J. Choi, D. Roundy, H. Sun, M. L. Cohen, and S. G. 
Louie, Phys. Rev. B 66, 020513(R) (2002). 

8 A. Y. Liu, I. I. Mazin, and J. Kortus, Phys. Rev. Lett. 87, 
087005 (2001). 

9 J. Kortus, 1. 1. Mazin, K. D. Belashchenko, V. P. Antropov, 
and L. L. Boyer, Phys. Rev. Lett. 86, 4656 (2001). 

J. W. Quilty, S. Lee, A. Yamamoto, and S. Tajima, Phys. 
Rev. Lett. 88, 087001 (2002). 

1 J. W. Quilty, S. Lee, S. Tajima, and A. Yamanaka, Phys. 
Rev. Lett. 90, 207006 (2003). 



K. P. Bohnen, R. Heid, and B. Renker, Phys. Rev. Lett. 
86, 5771 (2001). 

J.Hlinka, I. Gregora, J. Pokorny, A. Plecenik, P. Kus, 
L. Satrapinsky, and S. Benacka, Phys. Rev. B 64, 
140503(R) (2001). 

H. Martinho, C. Rettori, P. G. Pagliuso, A. A. Martin, 
N. O. Moreno, and J. L. Sarrao, Solid State Comm. 125, 
499 (2003). 

K. Kunc, I. Loa, K. Syassen, R. K. Kremer, and K. Ahn, 
J. Phys. Cond. Matt. 13, 9945 (2001). 
P. M. Rafailov, M. Dworzak, and C. Thomsen, Solid State 
Comm. 122, 455 (2002). 

D. DiCastro, E. Cappelluti, M. Lavagnini, A. Sacchetti, 
A. Palenzona, M. Putti, and P. Postorino, Phys. Rev. B 
74, 100505(R) (2006). 

T. Yildirim, O. Giilseren, J. W. Lynn, C. M. Brown, T. J. 
Udovic, Q. Huang, N. Rogado, K. A. Regan, M. A. Hay- 
ward, J. S. Slusky, et al., Phys. Rev. Lett. 87, 037001 
(2001). 

L. Boeri, G. B. Bachelet, E. Cappelluti, and L. Pietronero, 
Phys. Rev. B 65, 214501 (2002). 

J. Menendez and M. Cardona, Phys. Rev. B 29, 2051 
(1984). 



11 



21 A. Debernardi, S. Baroni, and E. Molinari, Phys. Rev. 
Lett. 75, 1819 (1995). 

22 G. Lang, K. Karch, M. Schmitt, P. Pavone, A. P. Mayer, 
R. K. Wehner, and D. Strauch, Phys. Rev. B 59, 6182 
(1999). 

23 M. Lazzeri, M. Calandra, and F. Mauri, Phys. Rev. B 68, 
220509(R) (2003). 

24 M. Calandra, M. Lazzeri, and F. Mauri, accepted Physica 
C (2007). 

25 A. Shukla, M. Calandra, M. d'Astuto, M. Lazzeri, 

F. Mauri, C. Bellin, M. Krisch, J. Karpinski, S. M. Kaza- 
kov, J. Jun, et al., Phys. Rev. Lett. 90, 095506 (2003). 

26 A. Q. R. Baron, H. Uchiyama, Y. Tanaka, S. Tsutsui, 
D. Ishikawa, S. Lee, R. Heid, K.-P. Bohnen, S. Tajima, 
and T. Ishikawa, Phys. Rev. Lett. 92, 197004 (2004). 

27 Y. Kong, O. V. Dolgov, O. Jepsen, and O. K. Andersen, 
Phys. Rev. B 64, 020501 (R) (2001). 

28 M. Calandra and F. Mauri, Phys. Rev. B 71, 064501 
(2005). 

29 E. Cappelluti, Phys. Rev. B 73, 140505(R) (2006). 

30 R. Verbeni, F. Sette, M. Krisch, U. Bergmann, B. Gorges, 
C. Halcoussis, K. Martel, C. Masciovecchio, J. Ribois, 

G. Ruocco, et al., J. Synchrotron Radiation 3, 62 (1996). 

31 C. Masciovecchio, U. Bergmann, M. Krisch, G. Ruocco, 
F. Sette, and R. Verbeni, Nucl. Instrum. Meth. B 111, 
181 (1996). 

32 C. Masciovecchio, U. Bergmann, M. Krisch, G. Ruocco, 
F. Sette, and R. Verbeni, Nucl. Instrum. Meth. B 117, 
339 (1996). 

33 A. Bergamin, G. Cavagnero, and G. Mana, J. Appl. Phys. 
82, 5396 (1997). 

34 F. James, MINUIT Function Minimization and Error 



Analysis, Computing and Networks Division, CERN, 
Geneva (1998), (CERN D506 Long Writeup). 

35 M. Cardona, in Light Scattering in Solids II, edited by 
M. Cardona and G. Giintherodt (Springer, Berlin, 1982), 
vol. 50 of Topics in Applied Physics, chap. 2, p. 19. 

36 X. K. Chen, M. J. Konstantinovic, J. C. Irwin, D. D. 
Lawrie, and J. P. Franck, Phys. Rev. Lett. 87, 157002 
(2001). 

37 S. Baroni, S. de Gironcoli, A. D. Corso, and P. Giannozzi, 
Rev. Mod. Phys 73, 515 (2001). 

38 Baroni et al., PWSCF/espresso code is available at 
http: / /www.pwscf .org 

39 J.P.Perdew, K.Burke, and M.Ernzerhof, Phys. Rev. Lett. 
77, 3865 (1996). 

40 N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 
(1991). 

41 S. G. Louie, S. Froyen, and M. L. Cohen, Phys. Rev. B 26, 
1738 (1982). 

42 S. de Gironcoli, Phys. Rev. B 51, 6773 (1995). 

43 P. B. Allen, Phys. Rev. B 6, 2577 (1972), P. B. Allen and 
R. Silberglitt, Phys. Rev. B 9, 4733 (1974). 

44 D. Waasmaier and A. Kirfel, Acta Cryst. A51, 416 (1995). 

45 A. Q. R. Baron, H. Uchiyama, R. Heid, K. P. Bohnen, 
Y. Tanaka, S. Tsutsui, D. Ishikawa, S. Lee, and S. Tajima, 
Phys. Rev. B 75, 020505(R) (2007). 

46 M. Lazzeri and F. Mauri, Phys. Rev. Lett. 97, 266407 
(2006). 

47 P. Allen, in Metals, superconductors, magnetic materials, 
liquids, edited by G. Horton and A. Maradudin (North- 
Holland Publishing Company, Amsterdam, New York, 
1980), vol. 3 of Dynamical Proiperties of Solids, chap. 2, 
p. 95. 



